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Abstract. We show that the pointer basis distinguished by collisional 
decoherence consists of exponentially localized, solitonic wave packets. Based 
on the orthogonal unraveling of the quantum master equation, we characterize 
their formation and dynamics, and we demonstrate that the statistical weights 
arising from an initial superposition state are given by the required projection. 
Since the spatial width of the pointer states can be obtained by accounting for 
the gas environment in a microscopically realistic fashion, one may thus calculate 
the coherence length of a strongly interacting gas. 
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1. Introduction 

If a quantum system is brought into contact with an environment, the states forming 
its Hilbert space are no longer equals among each other. Rather, a particular set of 
wave functions, the pointer states, is singled out by the environmental interaction pQ. 
These states are distinguished by the fact that they retain their purity for a relatively 
long time, say, under a Markovian master equation describing the environmental 
influence, while superpositions of different pointer states decay rapidly into a mixture. 
The associated decoherence rate can be many orders of magnitude faster than any 
other relevant time scale, including that of the dissipativc effects induced by the 
environmental coupling [H [3] . 

The concept of pointer states plays a fundamental role in explaining both the 
emergence of classicality and the operation of quantum measurement devices within 
the framework of quantum theory [H [2j [3]. In practice, knowing the pointer states 
and their time evolution allows one to directly specify the dynamics of an arbitrary 
initial state after the decoherence time, without solving the master equation. It 
is thus important to understand the emergence of pointer states in microscopically 
realistic environments, to characterize their form and dynamics, and to show that 
they constitute a basis of the system Hilbert space. 

Several strategies have been proposed to determine the pointer basis. Sorting all 
pure states according to their linear entropy production rate allows one to select the 
states with minimal loss of purity [J]. Another idea considers the solitonic solutions 
of a nonlinear equation related to the open quantum dynamics [SJ [7J [S]. Both 
approaches were shown to be largely equivalent if a particle is linearly coupled to a 
bath of harmonic oscillators [Jj . The pointer states are Gaussian wave packets in these 
linear models El [Jj [8] , and the complete, finite time decoherence in the Gaussian 
basis can be proved rigorously for a free particle [9] . 
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The linear coupling of the system coordinate to the bath assumed in these studies 
is often justified on the time scale of dissipation. However, it may misrepresent the 
decoherence process since it implies that the decoherence rate grows above all bounds 
for increasingly extended quantum states. This lack of saturation of the decoherence 
rate is unphysical for local interactions (an artifact of the unboundedness of the 
position operator); it predicts rates, e.g. for collisional decoherence measurements 
|10j . that are too large by many orders of magnitude. Moreover, previous studies did 
not demonstrate the conceptually important requirement that the statistical weight of 
an evolving pointer state is given by its initial overlap with the initial quantum state. 

In this paper, we study the emergence of pointer states due to collisional 
decoherence, allowing for a realistic and non-perturbative description of the impact of a 
gas environment [HIE]. We demonstrate how the orthogonal unraveling of the master 
equation [13j naturally explains the formation of (non- Gaussian) pointer states and 
their ensuing classical evolution, as well as the appearance of the expected statistical 
weights. Finally, by relating the widths of the pointer states to the microscopically 
defined localization rate of the master equation we can estimate the coherence length 
in an interacting thermal gas. 



2. The general pointer basis 

To define the notion of pointer states more precisely, consider a quantum system 
described by a Lindblad master equation dtp = C (p) due to its contact with an 
environment. One may say that C generates a pointer basis if the dynamics exhibits 
a separation of time scales, characterized by a decoherence time idee , such that there 
exists a unique set of pure states P a = |7r Q )(7r a |, which are independent of po and 
which have the property that at all later times any initial state po can be represented 
by a mixture of these states 

e £ Vo~ J daProb(a|po)Pa(t), if t > t dcc . (1) 

Crucially for interpreting the P a as proper pointer states, whenever po is a 
superposition of mutually orthogonal pointer states P^, the ensuing probability 
distribution Prob(a|po) = J2p W P ~ P) ^ s determined by the initial overlaps 

W/3 = Tr(poP/3(0)). (2) 

P a evolve in time, slowly compared to the decoherence scale idee, and they initially 
form a basis (usually overcomplete) with da being the associated measure, J da P a = I. 
The uniqueness of the evolving pointer states P a is a result of equation {TJ) applying 
to times t, which (while exceeding the decoherence time) can be much shorter than 
the above-mentioned time scale of dissipation, tdcc <C t <C idiss- 

Below, we demonstrate that a viable way of obtaining the P a and their time 
evolution is to follow the approach suggested in [5j El El [8] of using the solitonic 
solutions of a nonlinear dynamic equation in the space of pure states, 

d t P= [P,[P,£(P)]]. (3) 

This equation is satisfied by the pure state solutions of dtp — C(p), if there are 
any, while its nonlinear character is essential to distinguish pointer states from their 
superpositions. It has a natural interpretation as the deterministic part of a specific 
class of unravelings of the master equation, as explained below. 
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Figure 1. The nonlinear equation J7} drives any initial state |i/>o) (here a 
superposition of wave packets travelling to the right) into a localized soliton |7rt) 
that moves with fixed envelope and constant velocity. These solitonic solutions 
form an overcomplete set, the pointer basis of collisional decoherence. 



Having a solution P t to ([3]) at hand, one may generate further solutions if 
the system has an underlying symmetry. In particular, if C has the form C (p) = 
1/ (ih) [H, p] +D (p) and U( is a family of unitary operators satisfying 

U £>(/?) U f =C(UpU + ) , (4) 
ihd t U =[H,U], (5) 

then UfPtUj can be shown to solve ((3]). This will allow us to prove that the pointer 
states form a basis, rendering a single pointer state sufficient to construct the whole 
set. 

3. Collisional decoherence 

We now apply these ideas to the ID motion of a free particle subject to collisional 
decoherence [TTJ [121 HU • The corresponding master equation is the short time version 
of the full, non-perturbative, Markovian description of a test particle in a gas [15]. As 
such, it describes the loss of coherence but no dissipation effects, and it has the form 
C(p) = (2mih)^ 1 [p 2 ,p] + V(p), with incoherent part 

/oo 
dqG (q) e lq * /h pe-' lcl *l h - 1P . (6) 
-oo 

It involves the collision rate 7 and the normalized momentum transfer distribution 
G (q) , each defined by the scattering cross section and the gas temperature [Tl] [12]. 
and has Lindblad structure with bounded jump operators L q = (q)e lqx / h . If 

one disregarded the coherent part, an exponential decay of the position off-diagonal 
elements would ensue, (x\p t \x') = cxp[— F (x — x')t](x\po\x') , with a localization rate 
F (s) = 7 — 7 J dqG (q) exp (iqs/h) , which saturates at 7 for ,s — > 00. 

The form of C turns the vector representation of ([3J into a nonlinear integro- 
diffcrcntial equation 

dtipt (x) = - -^-.dl^t (x) - & (x) A [|Vt| 2 ] (x) , (7) 



dyg(y)g*F(y), (8) 

-OO 
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involving convolutions g * F (x) := J dy g (y) F (x — y). Since the dispersion caused 
by the first term in ([7]) competes with the localizing effect of the second term, all 
solutions of ([7]) turn into stable solitons, \ip t ) —> \^t) f° r t ^> t<jcc- These move with 
constant shape and velocity. |7r t+T (x)| = \wt (x — vqt) \ , as demonstrated numerically 
in figure [TJ Below, we show that these solitonic solutions can be identified with the 
pointer states of collisional decoherence. 

4. Soliton basis 

Even though the localization rate F(x) remains bounded for large x, the solitonic 
solutions to (J7J arc exponentially localized, n t (x) ~ e -k\x—v t\ e %tp(x—v t,t) as x -(- 00j 
with k > 0. This follows analytically by noting that for large x Eq. (J7J) takes the 
asymptotic form 

d t ipt (x) ~ - -^—dlipt (x) - (7 - cu) ip t (x) , as |x| -> 00 , (9) 
2mi 

with a constant = dy\tp t (y) \ 2 \tp t \ 2 * F (y) satisfying < < 7 . Here, we 

have used that the convolution \i/j\ 2 * G(x) vanishes for |x| — * 00, where G(x) denotes 
the Fourier transform of G(q). Any solitonic solution 7r t (x) = / (x — uoi) e lv ^ x ~ v ° t ^ 
thus takes the form / (x) = e~ fc ' x ' with fc > and a phase <p(x, t) asymptotically linear 
in x. 

Once a particular soliton P = 1 7r) (7r | has been found, further solitonic solutions 
to ([JJ are obtained by means of Galilei transformations. The phase space translation 
operator \J t = T Sj „ = exp (i (u t x — s t p) /H) satisfies condition ((4]) with ([6]). Moreover, 
condition ([5]) is met with H = p 2 / (2m) provided St = uot/m + so and ut = uq, 
rendering T SjU P T\ u also a solitonic solution of (J7J. This family of states, parametrized 
by r = (sq,uq), forms an overcomplcte basis since the identity, like any Hilbcrt- 
Schmidt operator, can be represented as I = J dTf(T)T s u QT^ u if Q is a trace-class 
operator [T^] . 

5. Orthogonal unraveling 

To identify the solitons as pointer states we now employ the method of quantum 
trajectories [TS1 U7J [TB] . In this framework, a pure initial state Po is propagated by 
a stochastic differential equation to generate an ensemble of pure states {PJ};, whose 
average yields the solution of the master equation, E[PJ] = e £t Po- Since there are 
infinitely many different stochastic processes with this property (unravelings), the 
ensemble corresponding to a fixed initial state has no observable consequences, apart 
form this average. 

However, Eq. ([1]) is a statement about the solutions for all initial states, such that 
the physically distinguished basis in |T]) may be related to a specific unraveling. In 
particular, a preference of a basis is obtained if the deterministic part of the unraveling 
exhibits solitonic solutions, such that the stochastic part no longer affects a trajectory 
once it has reached these solitons. 

We will now show, for the case of collisional decoherence, that this behavior is 
found with the orthogonal unraveling |13|j . Here, the wave function evolution has a 
deterministic part given by the nonlinear equation (|3|) , which is interrupted by random 

I Rather than the original version |13| . we use here the simpler orthogonal unraveling mentioned in 
|20l . It has the advantage that the jump operators are given explicitly. 
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jump^§|. The latter occur with the rate r q = (L*L 9 ) — (L*)(L g ) and are effected by the 
(nonlinear) operators J g = (L q — (L q )) / ^/rq~. 

First, we consider only the deterministic part of the quantum trajectory 
evolution, given by ([7]). We will restrict the discussion to pure initial states 
|"0o) = ^2iLi Ci\4>i), which are superpositions of narrow wave packets, with variance 
Var^x) < 2ttH 2 /V&raiq), situated on non-overlapping domains, fa (x) 4>*j^i (x) = , 
with separation Ax sufficiently large such that F (Ax) ~ F (oo) = 7. Under this 
assumption, which will be justified at the end of the following section, one can 
reformulate the nonlinear evolution |(7J) such that the essential part of the dynamics 
is already covered by the coefficients a (t), obeying a closed system of equations, 



iY 



-2 7 (^|c J | 4 -|c i | 2 ] N 2 . (10) 



df 



Id 2 



This equation was studied in |19| in the context of a model for state vector reduction. 
As shown there, all its stable fixed points have the form |cj| = S^ n , and a distinguished 
fixed point is approached monotonically, |c; (t — > 00 )| = 5i_ rn , with m = argmax^ ( | Ci | ) 
the index of the the largest wave packet contribution to \ipo}- 
The dynamics of the wave packets \fa) is described by 

d t fa (x) = - ^d 2 x fa (x) - fa (x) A [\fa\ 2 ] (x) 



N 



<(*) E 



\cjf% 0) , (11) 



with A defined in {8} and a rate 

%(x) = \fa\ 2 * F(x) - |0, | 2 * F(x) + 7 , (12) 

which is of the order of 7. In order to verify (p~0|) and (jlip . one takes the derivative of 
^(x) = ^2iCi(t)fa(t) with respect to time; using 

dy\fa(y)\ 2 \fa\ 2 *F(y)= 1 (l-5 tJ ) 7 (13) 

which follows from the above assumptions on fa(x), one finds that ipt(x) evolves 
according to 0. 

Since the coefficients Cj^ m tend to zero according to (jTUJ) , the coupling term 
vanishes after some time in Eq. (fTT|) . turning the latter into the soliton equation Q. 
In the absence of jumps, the initial superposition state thus evolves into the 
soliton which is associated to \fan), the wave packet contributing largest to |^o) (see 
Fig. 1). 

Next, we consider the probabilistic part of the quantum trajectories . Using the 
above assumption of localized wave packets, Var^ (x) < 2nh 2 /Varg (q) , one can 
approximate J dxe lqx ^ h \fa (x)\ 2 ~ e *<i x i/ h with xi = (x) ( ^ j . Under this assumption, 
the jump operator ] q affects only the coefficients. They are transformed as 

4 (?) = ^ U 9Xk/h £ e^A c k , (14) 



§ We note that Eq. ^ appears also in the context of a diffusive (i.e. not piecewise-deterministic) 
unraveling. It was used in [7] to study pointer states within a linear coupling model. 
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with Af being the normalization. The jump rate takes the form 




It vanishes for \a\ = 8i : „ , that is, when a stable fixed point of (TIT)]) is approached. The 
associated wave function can then evolve into a soliton state without being further 
perturbed by jumps. This confirms that the quantum trajectories of the orthogonal 
unraveling evolve into pointer states, the solitonic solutions to |Zj). 

Concerning the statistical weights, we now show how the stochastic process, 
described by (TTUj) . (TH| . (fTS)) can be treated analytically for N — 2 . The average effect 
of the jump (c' k )G = / dqG (q) c' k (q) simply interchanges the value of the coefficients, 
i.e., Kc'^gI 2 = | C2 1 2 - The probability of a trajectory starting from |ci (0)| 2 < 1/2 
to converge to |cj| = 8i t \ is therefore equal to the probability for an odd number of 
jumps. It is given by (1 — e~ 2 ^w)/2 f or this inhomogeneous Poisson process, with 
fi (t) = J Q drrtot (t) the integrated jump rate. In our case, 

d g r g = 2 7 | Cl | 2 (I-| Cl | 2 ), (16) 

where c\ (t) is a solution to l|10j) . Rewriting (fT0|) as 

2 7 | Cl | 2 (I-| Cl | 2 ) = i0 t ln(I-2| Cl | 2 ), (17) 

we find that fi (t — > oo) = — ln(l — 2 \c\ (0)| 2 )/2 . The probability for an odd number 
of jumps is therefore \c\ (0)| 2 = |(-0o|0i)| 2 - If \<Pi) = It,} this demonstrates that the 
statistical weights of the pointer states in the decohered mixture is given by the overlap 
K^oKi)! 2 , in accordance with 

The generalization of this result to N > 2 can be verified numerically. Taking 
G (q) to be a Gaussian with variance er^., we have implemented the stochastic process 
(jfOp . ([M]) . (Tl"5]) for various 2 ^ N ^ 100, using a Metropolis-Hastings algorithm to 
draw the momentum transfer q in accordance with the rate (|15|) . The initial states 
|-0o) were generated randomly by simplex picking. A x 2 -test then confirmed that the 
asymptotic trajectories, i.e., the pointer states P, = |7Tj)(7Tj|, arc distributed according 
to |<V>ok)| 2 ^TrlpoPi]. 

This suggests that the solitonic solutions to ([7]) are the pointer states of collisional 
decoherence in the sense that they retain their purity, while their superpositions 
decay into mixtures with weights given by the initial overlap. However, this result 
was derived under the assumption that the soliton variance is small compared to 
27rft 2 /VarG(q). We will see below that this assumption can be justified for small 
k = o^/mtvy. Still, the finite width leads to a small but finite loss of purity of the 
solitons, such that Eq. ((T|) is not an exact equality. 



6. Dynamics and size of the pointer states 

The motion of the pointer states can be characterized by their position and 
momentum expectation values also in presence of an external potential. Our numerical 
investigations indicate that the dynamics of this phase space trajectory changes from 
a quantum behavior to that expected from the corresponding classical mechanics for 
decreasing position spread ov of the pointer states. To explain this fact, we note that 
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an additional linear potential V (x) = ax in Q implies that the solitonic solutions 
have the form 7r t (x) = f (x — x t )e % ^ x ~ Xt ' t ' , with x t = —at 2 /2m + v t and / still 
exponentially localized. The pointer states thus move along the accelerated classical 
trajectories for linear potentials, implying that they will follow the general classical 
motion if the linearization of the potential is permissible over their spatial extension. 
Similar observations were made in 12 lj with linear coupling models. 

The pointer state position spread o~ n , which depends on the choice of G (q) and 
7, thus serves as an important quantity to characterize the mixture ([1]). We take 
G (q) to be a centered Gaussian with variance Oq. The dimensionless form of the 
soliton equation ([7]) then contains the ratio k = a^/mh^ as the only parameter. By 
numerically solving we extracted the dimensionless pointer width a^ac/h as a 
function of k. We found that the functional form can be reproduced by a simple 
localization model, which is inspired by studies of collisional decoherence [El HU H2] 
(and easily extended to the 3D case below): one assumes that scattering events 
occurring with rate 7 localize the wave function to a length scale £ loc , characterized 
by 1 — F(£i oc )/F(oo) =: exp (— a 2 oc /2), while it disperses freely between the collisions. 
Averaging the wave function width over the waiting time distribution of a Poissonian 
process then yields 

0-tt— = ha loc . (18) 

h 4aioc 

A value of a\ oc ~ 0.4 reproduces the numerically obtained widths very well (better than 
10%) over the full range of k. It follows from (fT8|) that for small k <C 4aj 2 oc , the soliton 
width is given by a 2 ~ a 2 oc H 2 /aQ, such that the assumption of small position variance 
of 4>i(x) can be justified for small k. In particular, we have checked numerically the 
approximation J dxe lqx ^ h \<fii {x)\ 2 ~ e %qXi ^, by using the solitonic solution of ([7|); the 
relative error is less than 2% for q G [— 2aa, 2ctg] and k ^ 10 -3 . 



7. Extension to 3D 



Equations (fTj)- (|15p are trivially extended to the 3D situation, though their numerical 
treatment is then more difficult. However, the above localization model allows one to 
directly estimate the 3D pointer width, and thus the coherence length of an interacting 
gas, using the microscopically realistic localization rate F derived in [T2]. We illustrate 
this for the case of hard-sphere s-wave scattering off a thermal gas, characterized by 
the thermal wave length A t h = y/2ir Jmk^Th. In this case, the localization scale 
60c = AocMth is determined by £ loc = exp (of oc /2 - 4tt^ 2 oc ) crfi (20r£i oc ) /4, which 
implies £i oc ~ 0.1 if we assume a\ oc = 0.4, as above. The average over the waiting 
time then yields the width of the pointer state 

ov = 7^7— + CiocAth, (19) 

where £f rcc is the mean free path. For a weakly interacting or thin gas, the pointer 
state width is thus essentially determined by £{ iee . In the limit of a strongly interacting 
or dense gas, on the other hand, it is bounded by the scale of the thermal wave length 
A th - 

We can use (U} to characterize the thermal state of a particle in presence 
of the interacting gas: the pointer state velocities display a Maxwell distribution, 
parametrized by A t h, while their position spread is given by (|19p . Identifying this 
state with the reduced single particle gas state /o gas , we can thus access the coherence 
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length A co h of the self-interacting gas. Since the latter is defined by the decay of the 
position off-diagonal elements, (x\p gaB \x') oc cxp(— tt \x — x'\ 2 /A^ oh ) , one obtains 

(20) 



with av from (|19[) involving the mean free path. This shows how the interactions 
in the gas reduce its ideal coherence length, a behavior possibly observable in the 
interference of cold, non-degenerate gas clouds [22]. 



8. Conclusion 



We have shown that the pointer basis of the collisional dccohcrence master equation 
is naturally obtained from the orthogonal unraveling, explaining the formation and 
dynamics of the pointer states, as well as the decay of an initial superposition into the 
correct mixture. Since collisional dccohcrence is a paradigm for the non-perturbative 
description of environmental influences, one may expect this unraveling to gcnerically 
provide the pointer basis. Apart from its conceptual importance, it also provides 
an efficient numerical treatment of open quantum dynamics, because the quantum 
trajectories cease to jump after the decoherence time, leaving the deterministic motion 
of the pointer states as the only remaining dynamics. 

We thank B. Vacchini for helpful discussions. The work was supported by the 
DFG Emmy Noether program. 
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